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Thorough control of quantum measurement is 
key to the development of quantum information 
technologies. Many measurements are destruc- 
tive, removing more information from the sys- 
tem than they obtain. Quantum non-demolition 
(QND) measurements allow repeated measure- 
ments that give the same eigenvalue . They could 
be used for several quantum information process- 
ing tasks such as error correction , preparation 
by measurement , and one-way quantum comput- 
ing . Achieving QND measurements of photons is 
especially challenging because the detector must 
be completely transparent to the photons while 
still acquiring information about them ' . Recent 
progress in manipulating microwave photons in 
superconducting circuits has increased demand 
for a QND detector which operates in the giga- 
hertz frequency range. Here we demonstrate a 
QND detection scheme which measures the num- 
ber of photons inside a high quality-factor mi- 
crowave cavity on a chip. This scheme maps a 
photon number onto a qubit state in a single- 
shot via qubit-photon logic gates. We verify the 
operation of the device by analyzing the aver- 
age correlations of repeated measurements, and 
show that it is 90% QND. It differs from previ- 
ously reported detectors ' because its sensitiv- 
ity is strongly selective to chosen photon number 
states. This scheme could be used to monitor the 
state of a photon-based memory in a quantum 
computer. 

Several teams have engineered detectors which are sen- 
sitive to single microwave photons by strongly coupling 
atoms (or artificial atoms) to high-Q cavities. This archi- 
tecture, known as cavity quantum electrodynamics (cav- 
ity QED), can be used in various ways to detect pho- 
tons. One destructive method measures quantum Rabi 
oscillations of an atom or qubit resonantly coupled to the 
cavity . The oscillation frequency is proportional to 
^/n^ where n is the number of photons in the cavity, so 
this method essentially measures the time-domain swap 
frequency. 

Another method uses a dispersive interaction to map 
the photon number in the cavity onto the phase differ- 
ence of a superposition of atomic states {\g) ^e^^\e)) / y/2. 



Each photon number n corresponds to a different phase 
^, so repeated Ramsey experiments can be used to es- 
timate the phase and extract n. This method is QND, 
because it does not exchange energy between the atom 
and photon. However, since the phase cannot be mea- 
sured in a single operation, it does not extract full in- 
formation about a particular Fock state \n) in a single 
interrogation. Nonetheless, using Rydberg atoms in cav- 
ity QED, remarkable experiments have shown quantum 
jumps of light and the collapse of the photon number by 
measurement. ' 

Here we report a new method which implements a set 
of programmable controlled-NOT (CNOT) operations 
between an n-photon Fock state and a qubit, asking the 
question "are there exactly n photons in the cavity?" A 
single interrogation consists of applying one such CNOT 
operation and reading-out the resulting qubit state. To 
do this we use a quasi-dispersive qubit-photon interac- 
tion which causes the qubit transition frequency to de- 
pend strongly on the number of photons in the cavity. 
Consequently, frequency control of a pulse implements a 
conditional tt rotation on the qubit - the qubit state is 
inverted if and only if there are n photons in the storage 
cavity. To ensure that this is QND, the qubit and stor- 
age cavity are adiabatically decoupled before performing 
a measurement of the qubit state. 

To realize this method we extend circuit-based cav- 
ity QED by coupling a single transmon qubit ' si- 
multaneously to two cavities. This allows one cavity to 
be optimized for fast readout and the other for coher- 
ent storage of photons. Related work by Leek et al. 
realized a single transmon coupled to two modes of a sin- 
gle cavity, where the the two modes were engineered to 
have very different quality factors. A schematic of the 
two-cavity device is shown in Fig. 1(a). A high-Q cav- 
ity serves as a photon memory for preparation and stor- 
age, and a low-Q cavity is used for fast readout of the 
qubit. The cavities are realized as Nb coplanar waveg- 
uide resonators with A/2 resonances at co's/27r = 5.07 GHz 
and co'rn/27r = 6.65 GHz, respectively. The cavities are 
engineered, by design of the capacitors Cg and Cm, to 
have very different decay rates {hZs/'^'Tr = 50 kHz and 
tZm/'^Tr = 20 MHz) so that the qubit state can be mea- 
sured several times per photon lifetime in the storage 
cavity. A transmon qubit is end-coupled to the two cavi- 
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FIG. 1. Circuit schematic and device, a, Circuit schematic 
showing two cavities coupled to a single transmon qubit. The 
measurement cavity is probed in reflection by sending mi- 
crowave signals through the weakly coupled port of a direc- 
tional coupler. A flux bias line allows for tuning of the qubit 
frequency on nanosecond timescales. b, Implementation on a 
chip, with cJrn/STT = 6.65 GHz measurement cavity on the left 
and its large coupling capacitor (red), and ujs/2n = 5.07GHz 
storage cavity on the right with a much smaller coupling ca- 
pacitor (blue). A transmon qubit (green) is strongly coupled 
to each cavity, with 5's/27r = 70 MHz and 5'm/27r = 83 MHz. 
It has a charging energy ^c/27r = 290 MHz and maximal 
Josephson energy Ej /2ti ^ 23 GHz. At large detunings from 
both cavities, the qubit coherence times are Ti ^ T2 ~ 0.7 /xs. 



ties, with finger capacitors controlling the individual cou- 
pling strengths (^s/27r = 70MHz and ^m/27r = 83 MHz). 
The usual shunt capacitor between the transmon islands 
is replaced with capacitors to the ground planes to re- 
duce direct coupling between the cavities. Additionally, 
a flux bias line allows fast control of the detunings 
A5 — ujg^e — and = ujg^e — between the trans- 
mon and cavities, where we use the convention of labeling 
the transmon states from lowest to highest energy as (^, 
e, /, h, ...). 

To achieve high photon number selectivity of the 
CNOT operations, there must be a large separation be- 
tween the number-dependent qubit transition frequen- 
cies. To obtain this, we use small detunings (A^/ Qs < 10) 
between the qubit and storage cavity. Figure 2 shows 
spectroscopy in this quasi-dispersive regime as a func- 
tion of flux bias when the storage cavity is populated 
with a coherent state ((n) ~ 1). Results of a numerical 
energy-level calculation are overlaid, showing the posi- 
tions of various transitions. We define cj^ g as the photon 
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FIG. 2. Pulsed spectroscopy with coherent state in storage 
cavity ((n) ^ 1) vs. qubit-cavity detuning As = ojg^e — 00 s. 
Calculated transition frequencies are overlaid in color. Red 
and orange lines are the \g) ^ \e) transitions of the qubit 
when n = and 1, respectively. Transitions to higher trans- 
mon levels (I/) and \h)) are visible because of the small de- 
tuning. The arrow indicates the flux bias current used during 
the CNOT operations. 

number-dependent transition frequency \n^g) 

Other transitions, such as |2,^) |0, /i), are allowed due 

to the small detuning. Fortunately, we also see that the 

separation between cj^ ^ and cj^ g grows rapidly to order 

^ 2g = 140 MHz as the qubit approaches the storage 

cavity. 

To test the photon meter, we generate single photons 
in the storage cavity with an adiabatic protocol. Our 
method uses the avoided crossing between the |0,e) and 
|1,^) levels to convert a qubit excitation into a photon. 
The preparation of a photon begins with the qubit de- 
tuned below the storage cavity (Ag — 3^s), where we 
apply a 7r-pulse to create the state |0,e). We then adi- 
abatically tune the qubit frequency through the avoided 
crossing with the storage cavity, leaving the system in the 
state 1 1,5'). The sweep rate is limited by Landau-Zener 
transitions which keep the system in |0, e). Our prepara- 
tion protocol changes the qubit frequency by 600 MHz in 
50 ns, giving a transition probability less than 0.1% (cal- 
culated with a multi- level numerical simulation). This 
protocol actually allows for the creation of arbitrary su- 
perpositions of 10,5) and \l,g) by changing the rota- 
tion angle of the initial pulse. For example, if we use 
a 7r/2-pulse, after the sweep we end up in the state 
( 1 0, 5) + e*"^ 1 1 , 5) ) / >/2, where (j) is determined by the rota- 
tion axis of the 7r/2-pulse. One could also use a resonant 
swap scheme, which has been successfully used to create 
Fock states up to \n = 15). The method used here has 
the advantage of being very robust to timing errors. 

After the photon is prepared, the qubit frequency is 
adjusted such that Ag/gs — 5. At this detuning, the sep- 
aration between cj^ g and cj^ g is ~ 65 MHz. In Fig. 3(a), 
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FIG. 3. Single photon preparation and CNOT selectivity, 
a, Pulsed spectroscopy vs. Rabi angle of preparation pulse, 
showing the reflected phase of a pulse at the measurement cav- 
ity frequency after a ^ 80 ns pulse near the qubit frequency. 
Traces are offset vertically for clarity and labeled with the ro- 
tation angle of the control pulse used in the preparation step. 
The dips correspond to c<;° e ~ 5.47 GHz and u;^ e ~ 5.41 GHz, 
respectively, b and c, Rabi driving the qubit transitions after 
preparing |n = 0) (b) and |n = 1) (c). The red (blue) traces 
show the measured qubit excited state probability after ap- 
plying an interrogation Rabi pulse with varying angle at cj^ e 
(cj^^e)- The residual oscillation of Ri{0) in c is mostly due to 
preparation infidelity. 



we show pulsed spectroscopy at this detuning for several 
rotation angles of the initial preparation pulse. We ob- 
serve well-resolved dips in the reflected phase of a pulsed 
signal sent at the measurement cavity frequency. The 
locations of these dips correspond to the qubit transi- 
tion frequencies for n = (co'^^g) and n = 1 (co'^^g), and 
the relative heights match expectations from the differ- 
ent preparation pulse rotations (e.g. a 7r/2-pulse results 
in equal height signals). 

To show selective driving of these transitions, we per- 
form Rabi experiments at cj^ ^ and cj^ ^ for the cases 
where we prepare |0,^) and |1,^). In each experiment 
we ensemble average measurements of the resulting qubit 
state after further decoupling the qubit from the storage 
cavity. For the \0,g) case [Fig. 3(b)] there is a large am- 
plitude oscillation when the drive is at cj^ ^ [red, Ro{0)] 
and almost no oscillation when the drive is at cj^ ^ [blue, 
Ri{0)]. When we prepare the situation is reversed 
[Fig. 3(c)]; however, in this case the residual oscillation 



of Ro{0) (red) is substantial due to small errors in the 
preparation of |1,^) associated with the initial rotation 
of the qubit and, more importantly, the ~ 10% proba- 
bility of energy decay during the subsequent adiabatic 
sweep through the cavity. 

The responses Ri{6) are a result of driving cj* ^ and 
the far off-resonant drive of co'^ g, where j ^ i. The cross- 
talk is seen in the small residual oscillation of Ri{0) in 
Fig. 3(b). In the supplement, we derive a method for ex- 
tracting a selectivity and preparation fidelity from these 
data, giving a selectivity > 95% for both interrogations 
and a preparation fidelity of \{n = llV^)!^ ^ 88%. These 
numbers were confirmed by doing equivalent experiments 
over a range of preparation pulse rotation angles between 
and 27r (not shown). 

If TT-pulses are used in the interrogation step, measure- 
ment results of the average qubit state directly correlate 
with the probability of being in the states \n = 0) or 
|n = 1). Details of the scaling needed to do this transfor- 
mation when the selectivity is < 100% are presented in 
the supplement. These are the desired CNOT operations 
of the photon meter. If we now insert a variable delay be- 
fore interrogating, we find that Pq (Pi), the probability of 
being in |n = 0) (|n = 1)) decays exponentially towards 1 
(0), as shown by the red (orange) trace in Fig. 4(b). The 
decay constant of Ti 3.11 ± 0.02 /is agrees with the 
linewidth of the storage cavity, 1/his = l/(27r50kHz) = 
3.18 /is, measured in a separate, low power (n ~ 1) re- 
flection experiment. 

Strong QND measurements are projective, such that if 
the measurement observable commutes with the Hamil- 
tonian, the system will remain in an eigenstate of both 
operators between measurements. Consequently, com- 
paring the results of successive interrogations provides a 
mechanism to test whether a particular protocol causes 
additional perturbations on the system. Here, we only 
compare ensemble average results, because the single- 
shot qubit readout fidelity for the device is ~ 55%. This 
is sufficient to reveal processes which change the photon 
number, and slight technical improvements to interro- 
gation speed or qubit readout fidelity should allow for 
real-time monitoring of the photon state. 

The protocol cannot be repeated immediately, though, 
because the first interrogation may leave the qubit in the 
excited state. To circumvent this problem, we use the fast 
decay rate of the measurement cavity to cause the qubit 
to spontaneously decay into the 50 Ohm environment. 
The "reset" protocol brings the qubit into resonance with 
the measurement cavity for a time, Treset = 50 ns, which 
is sufficient to reset the qubit with probability ~ 98%. 
The procedure is described in detail in Ref. 

After resetting the qubit, we can interrogate a sec- 
ond time. The full protocol for a repeated interroga- 
tion sequence is shown in Fig. 4(a). The combination of 
a CNOTo (CNOTi), a qubit measurement, and a qubit 
reset define an interrogation process Iq (/i). Data for 
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FIG. 4. Repeated measurements of photons, a, Experiment 
protocol. A microwave pulse and adiabatic sweep load a single 
photon into the storage cavity in the preparation step. This 
photon is interrogated repeatedly by number- selective CNOT 
gates on the qubit, followed by adiabatic decoupling, qubit 
readout, and reset, b. Single and repeated interrogation after 
preparing |n = 1) (top) or |n = 0) (bottom), ensemble aver- 
aged over ^ 50, 000 iterations. The near-perfect overlap be- 
tween single and repeated results demonstrate that the proto- 
col is highly QND. c. Transition probability diagrams for the 
interrogate n = (/o) and interrogate n = 1 (/i) processes. 
We extract 70 (71) = 1 (10) ± 3% and (^1) = 7 (3) =b 3%. 



the four possible combinations of interrogating |n = 0) 
and \n = 1) are shown in Fig. 4(b) as a function of delay 
between the first and second interrogations. The data 
are ensemble averaged over all results from the first in- 
terrogation, so we do not observe projection onto num- 
ber states. Instead, we again observe exponential decay, 
where the result of the second measurement is essentially 
indistinguishable from the first, indicating that the inter- 



rogation is highly QND. 

Deviations from the average measurements of a single 
interrogation stem from finite photon lifetime in the stor- 
age cavity and non-QND processes which cause transi- 
tions to other photon numbers [Fig. 4(c)]. Recording the 
second interrogation results for different delays allows us 
to subtract the effect of photon Ti and calculate the tran- 
sition probabilities for the Iq and Ii processes . In prin- 
ciple, /q and 1 1 can cause transitions to photon numbers 
outside of the n G {0, 1} manifold; however, the absence 
of statistically significant deviations from Pq -\- Pi = 1 
suggests that any such effects are negligible. Instead, we 
observe 70 (71) = 1 (10) ± 3% and ^0 (^1) = 7 (3) ± 3%, 
demonstrating that this protocol is highly QND. 

The protocol presented here is a fast and highly QND 
measurement of single photons, which we believe can be 
extended to detect higher photon numbers. It should be 
possible to demonstrate the projective nature of the inter- 
rogation and create highly non-classical states of light via 
post selection, and eventually with higher fidelity read- 
out it should be possible to observe quantum jumps of 
light in a circuit. 
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SUPPLEMENT 
Measured Voltage Scaling 

When the interrogation selectivity is less than 100%, 
we need to account for undesired rotations to correctly 
calculate the state probabilities from the measured volt- 
ages. The details of our calibration procedure follow. 

If we prepare |n = 0) or |n = 1) at time t = 0, when 
we interrogate at some later time there is an additional 
probability pd to decay, giving the density matrices 

P° = \9){9\®\0){0\, 

= \9){g\®{Pd\0){0\ + {l-Pdm{l\}, 

where indicates preparing state |i) at t = 0. We can 
model the interrogation pulses as operations which act 
on p*: 

Uo = Ry{n) |0)(0| + Ry{e') ® |1)(1|, 

Ul = Ry{e) ® |0)(0| + Ry{TT) ® |1)(1|, 

where e and e' are small angles. After interrogation, the 
integrated homodyne response is 
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where n G {0, 1} is the Fock state of the cavity, r G 
{0, 1, 1}, Vg {Ve) is the voltage measured when the qubit 
is in \g) (|e)), AV = Ve - Vg, and = \e){e\. 

By abusing the notation shghtly and treating e, e' as 
probabihties rather than rotation angles, we can calculate 
the PF" 



V, 



W\ = Vg 



9^ 



= Vg 



AV, 
AV-e, 

AV {e'{l-pd)+Pd), 
AV {{1 - pd) + epd) . 



We measure these five voltages in calibration experi- 
ments and invert the equations to find the parameters 
Vg^ AV, e, and e^ Note that this does not require 
perfect preparation fidelity because the model includes 
decay between preparation and interrogation pd which 
will also capture any fixed preparation infidelity. This 
gives the selectivities, {1 — e) and (1 — e'), as well as the 
preparation fidelity, (1 — Pd)- 

An unknown mixture of n = and n = 1 is character- 
ized by a single probability p. 



P=l5)(<?l®{p|0)(0| + (l-p)|l)(l|}, 
which produces the responses 

= Vg + AV{p + e'{l-p)), 
= Vg + AV{e-p+{l-p)). 

This leads to a simple rescaling to transform Wq and Wf 
into Pq and Pi 

^ W^' - (Vg + AVe') 
° AV{1 - e') 

W{ - {Vg + AFe) 
AV{\ - e) 



Pi 



Error Estimate 

The primary challenge in these experiments is obtain- 
ing sufficiently accurate and precise control of the qubit 
frequency to do high-fidelity operations. The narrow 
bandwidth pulses used in the CNOT operations means 
that even a few MHz error in frequency control results in 
a significant rotation error. We use deconvolution tech- 
niques similar to those described in the supplement of 
Ref. ; however, the fiux bias current response function 
drifts on a time scale of about one day, making it difficult 
to eliminate all classical control errors. Even after apply- 
ing corrections, there is a remaining spread of 2 — 3 MHz 
in the qubit frequencies over the various realizations of 



preparation to interrogation delay. This translates into 
a 2 — 3% error in the probability to find the qubit in 
|e) after applying a conditional 7r-pulse. The errors bars 
reported in the lower panel of Fig. 4(b) are due to this 
systematic error. 
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